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K~    SCATTERING  IN  DEUTERIUM  AT  100  to  250  Mev/c 
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Lawrence  Radiation  Laboratory- 
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Berkeley,    California 

May  25,     1959 

ABSTRACT 

An  investigation  of  K   -meson  elastic  and  inelastic  scattering 
at   100  to  250  Mev/c  is  reported.      Experimentally  determined  cross 
sections  were  compared  to  theoretical  cross  sections  from  the  two 
Dalitz  phase-shift  solutions  for  K    -p  scattering.      Within  statistical 
accuracy,    the  second  Dalitz  solution  is   found  to  be  correct. 
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INTRODUCTION 

For  several  years  the  Alvarez  group  has  been  engaged  in  the 
study  of  interactions  between  low-energy  K     mesons  and  protons.      In 
order  to  extend  this  investigation,    the  Lawrence  Radiation  Laboratory 
bubble  chamber  was  filled  with  deuterium  and  exposed  to  low-energy 
K     mesons  produced  in  the  Bevatron.      The  experiment  was  designed  to 

(a)  check  the  branching  rules  for  hyperon  production  deduced  from  the 

principle  of  charge  independence  in  strange-particle 

1 
interactions, 

(b)  investigate  the  scattering  of  low-energy  K     mesons  from 
deuterons. 

Dalitz  has  performed  a  comprehensive  analysis  of  the  low 

energy  K     meson    -proton  interactions  in  terms  of  a  short  range 

2 
S-wave  interaction    .        Because  of  statistical  limitations  in  the  experi- 
mental data,    there  are  ambiguities  in  the  solutions.      The  main  purpose 
of  this  investigation  is  to  try  to  use  the  scatterings  that  occur  in 
deuterium  in  order  to  determine  the  nature  of  the  correct    solution. 
By  use  of  a  simple  model  of  the  deuteron,    the  scattering  cross  sections 
may  be  estimated  from  a  knowledge  of  the  elementary  K     meson  — nucleon 
scattering  amplitudes.      We  are  led  to  two  distinct  predictions  for  the 
K     meson— deuterium  cross  sections,    one  of  which  should  agree  with 
the  experimental  data. 

As  a  by-product  of  this  investigation,    we  shall  be  able  to  perform 
a  separation  of  the  real  scattering  events  from  events  involving  hyperon 
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production tai|t  looking  superficially  like  scatterings*      This  separation 
is  necessary  so  that  a  meaningful  comparison  may  be  made  between 
the  experimental  hyperon  branching  ratios  and  those  deduced  from  the 
principle  of  charge  independence  in  strange-particle  interactions. 


-6. 
EXPERIMENTAL  APPARATUS 

A  450-Mev/c  K   -meson  beam  was  obtained  from  the  Bevatron 

through  the  use  of  a  target  and  spectrometer  system  reported  by 

3 
Horwitz  et  al.    (Fig      1).        After  emerging  from  the  coaxial  electro- 
magnetic spectrometer  the  K     mesons  were  slowed  down  by  a  suitable 
absorber  so  that  approximately  75%  came  to  rest  in  the   15-inch 
Lawrence  Radiation  Laboratory  deuterium  bubble  chamber  (Fig.    2). 
The  resulting  incident  beam  contained  one  K     meson  per   1000  back- 
ground tracks.. 

Four  stereo  views  of  each  expansion  were  taken  on  35  mm  film 
by  cameras  positioned  above  the  chamber.      After  a  scanning  and 
identification  of  events,    measurements  were  carried  out  with  the     UCRL 
precision  track-measuring  machine  (Franckenstein).      This  device 
-  semiautomatically  punches  onto    IBM    cards  the  pertinent  coordinates 
of  tracks  undergoing  measurement.      All  further  computations  are  then 
performed  with  the  aid  of  an    IBM  650  digital  computer. 

ANALYSIS  OF  DATA 
General  Analysis   Techniques 

Approximately  46,500  photographs  were  taken  during  the  run 
These  were   "pick-tooth"  scanned  by  trained  technicians.     In  this 
method,    the  scanner  is  limited,    by  means  of  a  small  mask  on  the  film 
retainer  of  the  scanning  projector,    to  examination  of  the  upstream 
10  cm  of  the  entering  beam       By  looking  only  for  incoming  negative 
heavily  ionizing  tracks,    the  scanner  obtains  a  sample  of  events  in- 
dependent of  the  nature  of  the  final  interaction.      Proper  normalization 
then  allows  calculation  of  branching  ratios  and  interaction  cross  sections 
free  of  scanning  biases. 

Of  the  3300  K~  mesons  that  entered  the  chamber,    2200  interacted 
to  produce  hyperons  and  450  decayed.      The  remainder  did  not  interact 
and  had  sufficient  momentum  to  pass  through  the  chamber.      Using  the 
scanners"    records,    a  physicist  verified  the  nature  of  each  interaction 
Useful  events  were  then  sketched  and  measured  on  the  Franckenstein. 
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Fig.    1.      Target  and  spectrometer  system. 
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Fig.    2.      Lawrence  Radiation  Laboratory  15-inch  deuterium 
bubble  chamber. 
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Selection  of  Events 

In  order  to  eliminate  biases  in  the  selection  of  events,    it  was 
necessary  to  limit  acceptable  events  to  a  volume  of  the  chamber,    the 
"fiducial  volume"  (defined  with  respect  to  the  fiducial  marks),    satisfying 
the  following  criteria: 

(a)  Since  a  good  measurement  of  the  momentum  of  the  incident 
K     meson  was  critical   the  distance  between  the  point  at 
which  the  meson  became  illuminated  and  the  point  at  which 
it  entered  the  fiducial  volume  must  be  greater  than  6  cm, 

(b)  Because  of  the  inhomogeneity  of  illumination  near  the  edge 
of  the  chamber,    the  fiducial  volume  must  be  confined  to  a 
distance  greater  than    5  cm  from  the  edges. 

(c)  The  fiducial  volume  must  be  greater  than  5  cm  from  the  top 
or  bottom  of  the  chamber  so  that  hyperons  produced  in 
interactions  decay  within  the  visible  volume  of  the  chamber. 

It  was  found  possible  to  select  a  satisfactory  volume  which  was 
still  large  enough  so  that  7  5%  of  all  interactions  occurred  within  it. 
The  lateral  extent  of  this  fiducial  volume  was  defined  by  a  template 
which  could  be  used  on  the  scanning  projector  so  that  events  lying  outside 
this  volume  could  be  immediately  rejected.      The  final  selection  of  events 
was  made  by  an    IBM    650  program  which  examined  the  coordinates  of 
the  interaction  vertex  with  respect  to  fiducial  marks  on  the  glass  window 
of  the  bubble  chamber. 

Kinematic  Analysis  of  Scatterings 

Each  event  was  treated  first  as  an  elastic  scattering.      For  the 
cases  in  which  the  K     meson  and  the  recoil  deuteron  came  to  rest  in 
the  chamber,   the  momentum  (magnitude  and  angle)  of  the  scattered  K 
meson  in  the  laboratory  system  were  taken  as  known  quantities.      An 
IBM  650  program  computed. 

(a)    the  momenta  of  the  incident  K     meson  and  the  recoil 
deuteron, 
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(b)  the  angle  between  the  incident     K     meson  and  the  recoil 

deuteron, 

(c)  the  angle  of  scattering  in  the  center-of-mass   system. 
Both  measured  and  computed  values  were  included  in  the  output  of 
this  program.      (See  Appendix  I.  ) 

For  a  second  analysis,    each  event  was  treated  as  an  inelastic 
scattering.      The  vector  momenta  of  the  scattered    K     meson  and  the 
recoil  particle  (now  assumed  a  proton)  were  used  as  input  data.      The 
program  computed 

(a)  the  momentum  of  the  incident  K~  , 

(b)  the  momentum  of  the  recoil  neutron, 

(c)  the  azimuthal  and  dip  angles  of  the  scattered  K     in  a  new 

coordinate  system, 

(d)  the  azimuthal  and  dip  angles     of  the  recoil  proton  and  neutron 

in  this  new  coordinate  system, 

(e)  the  azimuthal  and  dip  angles  of  the  neutron  in  the  laboratory 

system. 
(See  Appendix  II  for  discussion  of  problem.  ) 

The  momentum  transfer  to  the  deuteron  for  an  elastic  scattering 
was  computed  for  each  event. 

Classification  of  Events 

When  K     mesons  interact  in  deuterium,    many  possible  interactions 
may  occur.      Because  the  strangeness  quantum  number  associated  with 
K     mesons   (-1)  must  be  conserved,    all  absorptions  lead  to  the  production 
of  hyperons.      The  most  important  reactions  are: 


K 

-♦71        +71 

K" 

-*■    JJ-    +        V 

K" 

0 

-♦    (JL        +71        +      V , 
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K"  +  D  -  2+  +  tt"   +  n  (1) 

K"  +  D  -  2~   +  ti+  +  n  (2) 

K"  +  D  -  2°  +  7i°  +  n  (3) 

K"  +  D  -  2=  +  7i°  +  p  (4) 

K~  +  D  -  2°  +  tT   +  p  (5) 

K~  +  D  -  A      +  ti°  +  n  (6) 

K~  +  D  -  A     +  ti"   +  p            .  (7) 

In  addition,    the  K     meson  may  decay  in  flight, 

0 

(8) 

or  simply  scatter  from  a  deuteron  before  undergoing  an  interaction, 

K"  +  d  -  K"  +  d  (elastic),  (9) 

K~   +  p  +  n  (inelastic).  (10) 

We  are  concerned  with  the  reactions   (9)  and  (10)  in  this  work. 
Unfortunately,     many  of  the  reactions   (1)  through  (8)  have  the   same 
appearance  in  the  bubble  chamber,    and  frequently  it  is  possible  to 
assign  only  a  most  probable  interpretation  to  a  given  event. 

Before  a  detailed  kinematic  analysis  was  attempted  all  events 
relevant  to  this   study  were  separated  into  two  groups: 

1.  Definite  scatterings   (Fig,    3) 

(a)  All  interactions  followed  by  the  production  of  a 
definitely  identified  hyperon  (we  thus  know  that  a 
K     meson  left  the  first  vertex). 

(b)  All  interactions  for  which  the  length  of  the  emergent 
particle  was   greater  than  2.5  cm.      (This  ensured 
that  reactions  of  Type   (4)  were  not  included,    since 
it  is  extremely  improbable  that  the  low-momentum 
2     hyperons  produced  by  absorption  would  travel 
2.5  cm,    before  undergoing  decay), 

2.  Possible  scatterings  or  possible     2   n    p.      These  events  may 

be  subdivided  according  to  the  fate  of  the     2     hyperon. 
(a)     2     tt    p  events   (Fig,    4).      For  these  the    2     hyperon 
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ZN-2152 


Fig.    4.      Possible  scatter,    possible      2     tt  p  with     2     decay 
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decays  in  flight, 

2    -*  tt     +  n  . 

(b)  2  «    tt   p    events   (Fig.    5).      Here  the     2     comes  to  rest 
and  gets  captured  by  a  deuteron, 

2~   +  d    —2"   +  n  +  n 

—  A  +  n  +  n  , 
and  the    Ahyperon  decays  via  the  charged  mode, 
A   ->  tt"  +  p  . 

(c)  2      tt    p  events   (Fig     6).      These  events  are  similar 

r 

to  Type   (2)  except  that  the     A  decays  via  its 

neutral  mode, 

A  -tt°  +  n  . 
We  next  separated  the  definite  scatterings  into  elastic  and 
inelastic   scatterings.      If  we  examine  the  kinematics  plot  for  elastic 
scattering,    we  see  that  it  can  be  roughly  divided  into  three  regions  of 
recoil-deuteron  momentum  (Fig.    7).      For  elastic   scatterings  in  region 
A,    the  recoil  deuteron  is  not  visible  in  the  bubble  chamber  because  of 
its  low  momentum.      On  the  other  hand,    the  deuteron  is  always  visible 
for  elastic   scatterings  in  region  C.      If  the  elastic   scattering  falls  in 
region  B,    the  deuteron  may  or  may  not  be  visible,    depending  on  its 
dip  angle.      On  the  basis  of  a  comparison  between  measured  and  computed 
values  for  both  the  elastic  and  inelastic  hypotheses,    each  definite 
scattering  was   classified  as  elastic,    inelastic,    or  indeterminate.      The 
indeterminate  scatterings  were  further  broken  down  by  assigning  a 
possible  and  probable  classification.      The  criteria  used,    in  order  of 
relative  weight  assigned,    are: 

a.  momentum  and  (or)  length  of  recoil  deuteron, 

b.  angle  between  incident  K     and  deuteron, 

c.  coplanarity  of  event, 

d.  momentum  transfer  to  the  deuteron, 

e.  momentum  of  incident  K    . 

This  classification  of  events  is   shown  in  Table  I 
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Fig.    6.      Elastic  scatter  kinematics  plot  . 
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Table  I 


Definite- scattering  classification 


Elastic  112 

Inelastic  42 

Indeterminate  12 

a.  Probable  elastic        3 

b.  Probable  inelastic     3 

c.  Equally  probable        6 

Total  176 


The  determination  of  an  interaction  as  an  elastic  or  inelastic 
scattering  was  essentially  a  subjective  question  of  relative  goodness  of 
fit.   In  general,    whenever  an  event  could  be  fitted  to  the  elastic- scattering 
kinematics  within  limits,    it  was  assumed  to  be  elastic,    although  it 
obviously  always  fits  the  inelastic  kinematics  better.      Therefore  by 
relaxing  our  fitting  criteria  slightly,    we  could  combine  the  elastic  and 
probable  elastic   scatterings  in  a  single  group.      We  note  that  this  re- 
classification of  indeterminate  scatterings  into  elastic  and  inelastic 
scatterings  is  in  the  same  ratio  as  the  initial  elastic-to-inelastic 
classification. 

It  was  then  necessary  to  estimate  the  number  of     Sup  events 
in  the  possible  scattering,    possible       2   n   p     group.      Real  scatterings 
that  contribute  to  this  group  arise  from  the  scatterings  followed  by 
Reaction  (6)    or     K~  decay.      To  estimate,  the  contribution  for  Reaction  (6) 
we  may  extrapolate  from  the  number  of  observed  scatterings  followfed 
by  well-identified  events.      A  tabulation  of  all  scatterings  according  to 
length  of  the  scattered  track  and  final  interaction  is  shown  in  Table  II. 
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Table  II 


Length  of  scattered  track 


Event  Ending  L,       2.5  cm.  L-,       2.5  cm 


A  t.°  n                                                                           23  15 

0     a 

(A)  n     n                                                                            10  9 

K"  decay                                                                         13  40 

K~  through                                                                    21  0 

A  T."p                                                                             17  7 

(A)  v~p                                                                        9  5 

2       Ti"p                                                                      13  14 

1                                                     .-  2 

0  1 

6  6 

5  1 

4  2 

13  3 


""(A)     indicates  the  decay  was  the  neutral  mode 

b 
event  following  scattering  in  question  was  a  possible  scattering. 
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We  select  the    Att~p,    2     tt   n,    2   tt   n,    and  2   tt   n  interactions.      Of  the 
956  events  of  these  types  observed  in  the  fiducial  volume  of  the  chamber, 
29  scattered  within  2.5  cm.    of  the  final  interaction.      Therefore,    for  the 
470   An   n  and      (A)tt   n    interactions  observed,    we  expect 


x  470  (    -^—   )  a  16±4 
956 

of  these  events  to  have  scatterings  within    2.5  cm.    of  the  interaction. 

This  leaves   24  -16  =  8      2   tt  p    events  in  this  group.      From  the 

A  -decay   branching  ratio,    we  expect     5  2~  tt  p  and  3  2  tt  p  interactions,  If  we 

make  the  same  calculation,    accepting  all  of  the   1700  observed  interactions 

except  the  Att   n,    (A)tt   n    and    K      decay  events,    we  obtain  from  Table  II 


*470  ( li )  «  16 

1700-470 

as  the  number  of      Att   n     events  with  scatterings   within  2.5  c.  m.    of 
the  interaction. 

We  also  expect  a  contribution  to  the  possible      2  tt   p  classification 
from  those  scatterings  which  are  followed  by  a  K     decay.      By  measure- 
ment,   we  find  that  the  average  path  length  in  the  fiducial  volume  of  the 
chamber  for  a  K     meson  before  undergoing  decay  was   8.4  c  m.      Since 
13  K     mesons  scattered  at  distances  greater  than  2.5  cm.    from  the 
decay  vertex,    we  expect 

— (13)  a  5 


8.4  -   2.5 


actual  scatters  within  2.5  cm.    of  the  vertex.      This  leaves  40-5=  35 
real  events  of  the  Type   (4).      The  predicted       2    tt    p  separation  is   shown 
in  Table  III. 
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Table  III 


Predicted     S~ n   p    separation 


Sup  Scatterings 


Possible     2   -n   p  35  5 

Possible     S     tt   p  5  10 


A 

» — 


Possible     Z    tt   p  3  6 


In  addition,    an  independent  decision  was  made  as  to  classification 
of  each  possible  scattering  —  possible      2   n   p     event  from  examination 
of  the  kinematic  fit  to  the  elastic-   and  inelastic- scattering  hypotheses. 
Some  of  these  interactions  did  not  fit  a  scattering  with  the     K 
momentum  available  in  the  incident  beam  and  could  be  immediately- 
classified  as     Sup     events.      In  general,    possible     Sup  interactions 
for  which  the  incident     K       meson     appeared  to  come  to  rest  were 
classified  as  definite     S    tt    p.      The  remaining  interactions  in  this 
group  were  considered  indeterminate  and  assigned     a  possible  or 
probable  classification  as  before.      Table  IV  summarizes  the  experi- 
mental    Sup  separation. 
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Table  IV 


Experimental   i   ti    o  seoaration 


Definite   2    tt   p 


Definite   2    r.    p 
P      * 
Indeterininate 


-    0 

a.  Probable   2    t,    n  1 

P      A 

b.  Probable  scatter  4 

c.  Equally  probable  2 


29 


Indeterminate  1  1 

a.      Probable     Z   tt   p  4 

'  b.      Probable  scatter  3 

c.      Equally  probable  4 

Definite    2a  t;    p  4 

Indeterminate  1  1 

0 

a.  Probable     Z«t    p  2 

b.  Probable  scatter  3 

c.  Equally  Probable  6 
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Once  again  the  fitting  criteria  were  relaxed  so  that  probable 
2   v   p    became  definite     Sup,    and  scatterings   (both  probable  and 
equally  probable)  became  definite  elastic  or  inelastic  scatterings  to  be 
separated  as  in  the  definite-scattering  group,      We  conclude  that  we  have 
35  2   it   p,    6     2  .  rr   p,    and    3   S   tt   p. 

Path -Length  Determination 

Since  only  the   scatterings  occurring  in  the  fiducial  volume   were 
accepted  for  analysis,    it  was   necessary  to  calculate  the  K     path  length 
in  this  volume  for  various  momentum  intervals   so  that  cross   sections 
might  be  computed.      An    IBM    650  program  was  devised  to  compute 

a,  the  momentum  of  each  incident     K     meson  at  its  entrance 

to  the  fiducial  volume  (P   ,         n        ), 

chamber 

b.  its  total  path  length  in  the  fiducial  volume. 

To  circumvent  the  problem  of  a  possible  scanning  bias  against  observing 
higher --momentum  K     mesons  that  did  not  interact,    only  tracks  entering 
the  fiducial  volume  with  momentum',  less  than  225  Mev/c  were  retained. 
This  is  the  minimum  momentum  for  a  K     meson  just  passing  completely 
through  the  fiducial  volume  without  interacting.      An  examination  of  the 
momentum  distribution  at  the  entrance  to  the  fiducial  volume  for  those 
scatterings  followed    -by  K     mesons  undergoing  decay  showed  two  that 
entered  with  momentum  less  than  225  Mev/c.      Since  the  path-length 
contribution  to  the  momentum  intervals  is  not  the   same  for  a  typical 
K     decay  as  for  the    other  K     interactions,    it  was  decided  that  this 
group  could  be  eliminated  without  seriously  affecting  our  over-all 

statistical  accuracy.      This  also  eliminated  from  the  data  the 

v-    0 

2.._  tt   p  events. 

Another  IBM    program  was  used  to  sum  up  the  K     path  length 

in  25-Mev/c  intervals  for       0  <C  P   ,         ,        <^  ZZb  Mev/c.      The 
'  — *-    chamber  —  ' 

Att    n  and    An    p     groups  were  taken  as  typical,    completely  measured 
groups  for  input  to  this  program.      The  summation  was  made  first  by 
using  momentum  via  curvature  or  momentum  via  range,    depending  on 
which  was  higher.      This   gave  a  minimum  path  length,    especially  at  low 
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momenta,    because  when  the  momentum  via  curvature  is  only  slightly 
higher  than  momentum  via  range,    the  path-length  contribution  to 
intervals  below  100  Mev/c    is  lost.      Then  the  summation  was  made 
using    momentum  via  range.      .This  gave  the  maximum  path  length 
because  it  treats  all  tracks  as  stopping,    and  in  reality  there  is  a  20% 
in-flight  contamination.      However,    it  is  considered  that  the  latter 
determination  is  more  nearly  correct,    and  the  accepted  values  were 
taken  between  the  two  calculated  path  lengths  but  nearer  the  momentum- 
via-range  value.      These  were  normalized  to  the  total  number  of  incident 
K     mesons  and  the  final  path  lengths  are  shown  in  Table  V. 

Table  V 

K    path  lengths 

P   ,  ,         interval  Path  length 

chamber  ° 


50  to  74  Mev/c  432  ±  102  c.  m. 

75  to  99  1050  ±  199 

100  to  124  2025  ±  224 

125  to  149  3205  ±  159 

150  to   174  4305  ±     49 

175  to  199  4310  ±  116 

200  to  224  2260  ±     80 
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THEORETICAL  CONSIDERATIONS 


For  elastic  scattering  from  a  single  particle  at  the  origin,    the 

4 
wave  function  at  any  point  is 


+(5)='e*p(i*r.  PJ  +  fce^ySEEiiS), 


where  the  first  term  is  the  incident  wave  and  the  second  term  is  the 
scattered  wave.      All  quantities  are  evaluated  in  the  center-of-mas  s 
coordinate  system.      If  we  consider  only     S-wave  scattering,    the 
scattering  amplitude  is 


£{6,  c|>)  =  T   (a  constant). 
If  the  scattering  center  is  not  at  the  origin,    the   scattered  wave  is 
given  by 


exp(iK 


r   -   r 


0 


) 


exp  (iK  •    rQ    ) 


where     r        is  the  location  of  the  scattering  center  and    exp(iK  ■    r_) 
represents  the  amplitude  of  the  incident  wave  at  the  scattering  center. 
When  considering  elastic  scattering  from  the  deuteron  (i.  e.  ,    the  deuteron 
does  not  break  up  during  the  collision)  we  may  obtain  a  first  approxi*rnation 
to  the  proton- scattered  amplitude  by  averaging  over  the  location  of  the 
scattering  center  of  the  deuteron,    as   shown  in  Fig.    8. 


roton 


deutron 
center  of 
mass 


neutron 


Wr) 


Fig.    8.      Geometry  of  scatterinj 
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In  this  approximation,   the    total  scattered  wave  from  the  proton  (located 
R 


at 


~2~ 


)  is 

exp  (iK 

T       /  exp  (iK  •   ?J 
p     I  U 


R 


Now  we  see 


R 

r  "    T" 


+D(R) 


dV 


R 


R 


,    2      -         it  R*2    ,     2 

Mr      -  r    •     —   +    -^     ) 


1   - 


r   •    R 


2  r' 


~  r 


r   ■    R 

2r 


so  that  the  scattered  wave  becomes 


exp  (iK  •    rQ  ) 


,xp  (iK  [r  -  4_5.]) 


VR> 


dVR, 


r  •    R 

where      is   small  enough  to  be  ignored  in  the  denominator  but 


not  in  the  exponential,    where  it  acts  like  a  phase  shift.      Since  we  have 
we  have  also 


1*    - 


T     exp  (iKr) 
P"        r 

„       exp(iKr) 
P  r 


;xp  (iK  •      — j-   )  exp  (-iK  ■     — j-  ) 


*D(R) 


dV 


R 


exp  (iq  •    -j-  )       ^D(R) 


dV 


R 


where     q  =  k  -  ky       and      k  is  the  incident     K       momentum  in  the     K    -d 
center-of-mass   system  and    k'   is  the  outgoing     K     momentum  in  the 
same  coordinate  system.      By  the  same  argument,    the  scattered  wave 
from  the  neutron  is 
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T       exp  (iKr)  / 


n  r 


) 


exp  (iq  •     -^-    ) 


*D(R) 


dV 


R 


and  the  total  scattered  wave  is 


,rr,  ™    v     exP  (iKr) 

(T     +  T    )    — — - 

n  p  r 


exp  (iq*.      —  ) 


*DW 


dV 


R 


The  integral  becomes 


R 


F(q)  =        exp  (iq    ^-         cos    0     ) 


*DW 


2    R2    dR    d^ 


=  2     tt    /  exp  (iq      —     cos    9) 


VR> 


2    RZdRd  (cos  6) 


and  we  see  that  we  have 


R 


:xp  (iq     -y-      cos  6    )     d     (cos  6  )  = 


xp  (iq   — j-   cos  0) 
R 


cos  0=1 
COS0~~1 


exp  (iq'  -j-   )      -  exp  (-iq      — ^—  2  sin  (q   —  ) 


R 

iq      — 


R 

1    T 


Thus  the  scattered  amplitude  is 


j    exp  (iKr) 

scat  n  p '  r 


|  4    T!Sin  (q   -^  ) 
R 

q  -T- 


*D(R) 


2  .R2  dR. 


In  general, 


da  2  ,  * 


■dTT=  r    *  ~- > 


scat    scat 


T     +  T 
n  p 


2       F(q) 

iTToT 
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It  is  now  necessary  to  evaluate  the  elementary  scattering 
amplitudes  in  the  K    -nucleon  system  in  terms  of  experimentally 
determined  phase  shifts.      For  S-wave  scattering  we  have 


ik    .    ,  =  const 
total 


sin  (Kr  +  6) 
Kr" 


=  const 


sin  Kr  cos   6  +  cos  Kr  sin  6 
Kr 


where     6     is  the  S-wave  phase  shift.      But  in  general  we  had 

,                 rjjr         \    ,    rp,  exp  (iKr) 
Y  =  exp  (iK  •    r)  +  T  — £— ■ , 

and  if  we  now  break  this  up  into  angular -momentum  states  and  take 
Only  the     S-wave  (  i   -  0)  part,    we  obtain 


J 


exp  (iK  •    r)  =  2     (2£   +  1)  i   j^Kr)  p^cos    Q  ) 


Therefore 


sin  Kr 
Kr 


i        sin  Kr       ,    .^^     cos   Kr  +  i  sin  Kr 
qj  =  +  KT    


Kr 


Kr 


1   +  iKT  ,.  KT  cos  Kr 

sin  Kr   +    ■ 


Kr 


Kr 


1+iKT 
Kr 


sin  Kr  +(- 


1  +  i  KT 


)  (KT    cos     Kr) 


But  we  also  had 


ljj  =  const 


cos   6 
Kr 


sin  Kr  +  tan  6  cos  Kr 


2  8. 


Identifying  terms,    we  get 


also 


tan  6  = 


KT 


1   +  i  KT 


tan  6  +  i     KT  tan  6  =     KT 


T  = 


tan  6 


K(l-i    tan  6)   ' 


(1   +  i  KT)  tan  6  =  KT, 

(1  +  i  KT)  sin  6  =  KT  cos  6, 

(1   +  i  KT)      exp(i6)-exp(-i6)l   =  iKT    [exp(i6)  +  exp(-i6)~| 

(1   +  i  KT)     exp(2i6)-f]  =  iKT  [exp(2i6)  +  fl 

exp(2i6)  -   (1   +  iKT)  =  iKT 

2  iKT  =  exp(2i&)  -   1 


T  = 


exp(2i6)  -    1 
2iK 


This  then  is  the  formula  for  the  scattering  amplitude,      T,      in  terms  of 
of  the  S-wave  phase  shift,      6. 

Actually  the  K    -p     system  is   a  combination  of  equal  parts  of 
the  two  isotopic   spin  states,    1  =  0  and  1=1, 


K 


1 


si   2 


[x  (i)  +  x  ioj\m 


In  the  zero-effective-range  approximation,    the  phase   shifts  are  related 
to  the  S-wave  scattering  lengths,      A_     and     A.,    by  the  formulas 


tan  6~   =  K  A„, 
tan  6,    =  K  A.  , 


where     Aq     and     A.     are  complex  because  of  the  absorptive  processes 
(hyperon  production,    etc.  )  that  can  occur.      From  the  experimental 

K    -p    data,    Daltiz  has  computed  two  possible  solutions  for  the  scattering 

5 
lengths  as   shown  in  Table  VI. 
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Table  VI 


AQ  (X  1 0 l  3  cm)  Aj  (X  1 0 l  3  cm) 


Soln  a  ±  0.20  +  0.76i  ±  1.62  +  0.38i 

Soln  b  ±  1.88  +  0.82i  ±  0.40  +  0.41i 


The  real  parts  of  both  solutions  have  the  same  relative  sign  but  the 
absolute  sign  is  not  determined. 

In  the  K   -neutron  system,    only  the    1=1   state  is  involved,    thus 

K"N  =  X    (1)  • 
When  we  relate  the  scattering  lengths  to  the  scattering  amplitudes 
through  . 

T     =     - 


1   +  i  KA 
we  can  find  the  scattering  amplitudes  for  deuterium.      Since  we  have 

T  T'+T0 

V  2  ' 


we  can  write 


Ti 


3  1 

T+T=—       T,+     —        T 

n  P  2  l  2  0 
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RESULTS 

Using  the  ideas  presented  in  the  preceding  section,    Gourdin  and 
Martin  have  calculated  the  theoretical  differential  cross   section  for 
elastic   scattering  in  deuterium.        In  the  laboratory  system  it  is 


da 


d(cos  6) 


=   2  ti   T      +  T 


F(q) 


F(0) 


d(cos  9C    m    ) 


d(cos     0) 


They  have  evaluated  the  effects  of  the  form  factor, 
the  equation 


F(q) 


F(0) 


,    using 


F(q)  = 


qr 


f 


0 


+D<r) 


sin 


-± r      dr 


Assuming  that  the  deuteron  is  adequately  described  by  the  Hulthen 
wave  function,    we  have 


1 


^D^  ;     ~F~        exp  (~PT)  "  exp  ("ar^| 


Gourdin  and  Martin  find 


F(q)  1         4ap(a  +  B)  -1      l6q   (g-a)    (g  +  a) 

F(0)  q  (p   -   a)  q    +4q      T/R-nA^^?  m+n  \^lu 


where     a  =  45.7  and  (3  =  7a.      A  plot  of 


F(q) 
F(0) 


(P-d)    +2(p+a)    !;+64ap(a+P) 


vs     q    is  presented  in 


T' 


Fig.    9-      This  general  curve  can  be  used  to  obtain  the  differential 
cross   section  at  particular  laboratory- system  momenta.      For    175 
Mev/c  (lab)  we  have  computed 
solutions   we  find 

2 


T      +  T 
n  p 


For  the  two  Dalitz 


Solution  a: 
Solution  b: 


T      +  T 
n  p 

T      +  T 
n  p 


=   3.45, 


=  2.95 
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The  resulting  differential  cross  section  for  elastic  scattering     is 
shown  in  Fig.     10. 

To  obtain  the  total  elastic  cross  section  we  have  evaluated 
the  integrals, 

/  ^>dn=2*    l|g       d  cos    0 


2  71 


F(q) 


F(0) 


T     +  T 
n  p 


Simpson'  s  rule  was  used  to  obtain 


F(q) 


F(0) 


d  cos 


d  cos  6    from  Fig.    11 


for   several  values  of  incident     K    -meson  laboratory- system  momenta, 
and  total  elastic  cross  sections  were  computed. 

A  plot  of  the  angular  distribution  of  the  experimentally  determined 
scatterings   shows  a  large  peak  in  the  forward  direction.      This  indicates 
that  the  Coulomb  scattering  effect  is  large,    and  since  it  is  not  easily 
calculable,    we  have  discarded  all  scatterings  with  cos  9-^   0.95      (lab) 
and  have  corrected  the  total  cross   sections  for  this  effect. 

In  an  analogous  manner,  using  the  equations  of  Gourdin  and 
Martin,  we  have  computed  the  differential  and  total  inelastic  cross 
sections . 

Because  of  statistical  limitations  in  our  data,    we  have  used 
path  lengths  to  predict  the  number  of  events  we   should  expect  to  see  in 
each  momentum  interval  for  Solution  a    or  for  Solution  b.      Thus  we 
have 


N  =  I   a   NQ  =  1   a   (.0642)  (6.03X  1Q23)    (10"2?) 


=    (3.85  X10      )  i   a   , 


where     <j     is  in    millibarns.      A  summary  of  the  predicted  values  is 
shown  in  Table  VII. 
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Fig.    9.      Form  factor  vs  cos  6 
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Fig.    10.      Total  scattering  cross-sections   - 
Solution  a. 


-33 


cvi_ 

LL 


1.0 

0.9 

0.8 

|    0.7 

Zl   0.6 

"     0.5 

04 

0.3 

0.2 

0.1 


1          1          1          1          1 

PK  =  25    Mev/C 

1                     1                     1                     1 

PK=5Q - " " 

-— — — " '^^^^^y^' 

__ " 

^^^~*^     ^ ~^r    y ///    - 

PK  =  75               "-" "                                                -_--— -- 

= "                                                          —■ — -^"■"" 

^~-"~**^'  ^^'^^S^s// 

PK=I00^ ^ 

^^^^\^^^jy/y// 

P K  =  125 '      "                                  '        ^^x*' 

^ ^      V^y// 

P  K  =  1  5  0__— - ' ""                   ^                 ^-^-"-"' 

'^\^^Z^ 

PK=I75 . """"""            __— - ""'    ^«— ■ \- 

„ ^^ 

"      pkT2£0_— — — " """H~— — ""-^^Z- — ^-"' 

-      p,T?25_ -—-^__- * 

—1 T 

1                      1                     1                    1                     1 

■1.0      -0.8       -0.6        -0.4 


0.2  0.0  0.2  0.4 

Cos    0|2  (cm  ) 


0.6 


0.8 


1.0 


MU-17502 


Fig.     11.      Total  scattering  cross  -  sections 
Solution  b. 
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Table  VII 


Predicted  cross-section  data 


P,    ,  interval  cr    ,  0"  •       , 

lab  el  inel 

(Mev/c)         li  d  a        b 


t 
4mb)_ 


Corrected 

<7t  Nt 


100  -   149  359       150  74      38 

150  -   199  210         92  73      39 

200  -   224  139         64  64      35 


433  188 
281  130 
202         99 


406    179       86  38 

262    120       91  42 

187      87       1_7  _8_ 

194  88 


We  can  compare  the  experimentally  determined  number  of 
scatterings  and  the  predicted  numbers  to  determine  which  solution  is 
correct.      If  we  accomplish  this  through  the  total  scattering  cross   section, 
our  elastic-   and  inelastic- scattering  separation  is  not  involved.      We 
find  58  scatterings  occurring  between  100  and  225  Mev/c  that  fit  our 
acceptance  criteria,    and  therefore  conclude  that  Solution  b  is  correct. 

In  summary,    we  present  the  elastic,    inelastic,    and  tota] 
scattering  cross   sections  for  the  two  solutions  in  Figs.     12  and   13.      We 
have  also  plotted  the  experimentally  determined  cross   sections  for 
P,    ,    =    125,     175,    and  212.5  Mev/c.      In  Fig.     14  we  have  normalized  the 
theoretical  curves  to  the  observed  numbers  of  events  in  order  to  compare 
the   shapes  of  the  angular  distributions. 


35 


500 


400 


£      300- 


200 


100 


oslQb<0.95) 


100  150  200  250  300 


b       (  Mev  /c  ) 
lab 


MU-17504 


Fig.     12.      Normalized  angular  distribution  at 
175  Mev/c. 
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Fig.    13.      Total  scattering  cross-sections, 
solution  6 . 
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Fig.    14.      Normalized  angular  distribution  at 
175  Mev/c 
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CONCLUSIONS 

Because  of  the  enormous  difference  in  the  predicted  K   -D 
cross  sections  for  the  two  Dalitz  phase-shift  solutions,    it  has  been 
possible  to  make  a  strong  argument  for  Solution  b,    even  with  the 
limited  statistical  data  available.      It  must  be  remembered  that  this 
conclusion  is  independent  of  our  ability  to  separate  elastic  and  inelastic 
scatterings,    since  only  the  total  cross  sections  were  used. 

In  addition,    comparison  of  the  elastic  and  inelastic  components 
of  the  scattering  events  with  the  theoretical  curves  shows  that  the 
simple  model  used  by  Gourdin  and  Martin  to  calculate  the  scattering 
in  deuterium  is  adequate  to  describe  the  angular  distributions,    and, 
approximately,    the  ratio  of  inelastic  to  elastic  events.      It  will  be 
noticed  however,    that  the  absolute  values  measured  are  consistently 
lower  than  those  predicted.      It  is  expected  that  more  detailed  theoretical 
calculations  which  include  the  effects  of  secondary     scattering  will 
result  in  lower  cross   sections  than  presented  here.      If  future 
calculations   show  that  the  scattering  amplitude  for  Solution  b     can  be 
depressed  by   '50%,    then  the  predicted  cross   sections  will  decrease 
by   100%  and  no  conclusion  can  be  drawn  from  this   experiment. 
Experience  with  the  effects  of  higher-order   (  multiple- scattering)  effects 
in  the       n-meson  nuclear   system  leads  one  to  expect  that  such  a  large 
depression  will  not  occur. 
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APPENDIX 


A.      Elastic  Scattering  Analysis  for  K~  +  d    -»K     -f  d 


Given:     P*2    and     Q 


1  =      incident     K 

2  =      scatter     K 

3  =     recoil    d 


Find:     P^ ,    P^,    0^,    0    2(c.  m.  ),    LC3(proj) 

Let     P.     be  in  the    x    direction  and  the  plane  of    P-,    and    P,     be  th< 
xy    plane.     From  momentum  and  energy  conservation  we  get 


P,    =  P-,  cos  0    ?  +  P_  cos  6,-.   , 


0      =  P.  si 


in 


12 


+  P0   sin  0 


13  ' 


1 


P    2  *S 

_1L_  +  JL 


(Al) 
(A2) 

(A3) 


§22 
mR  mK  md 

Solving  (Al)  and  (A2)  for     P3 

P3Z  =  P    2  -   2  PjP2  cos 

Multiplying   (A3)  by  2m 

p    2  2         ^K_  2 

1  2     +   m,  3 


+  P    ' 
12  +  *2 


2  2 

P2     +RP3 


(A4) 


(A5) 


Substituting  (A4)  in  (A5)  and  dividing  by     P2    ,    we  get  the  quadratic  eqn, 


P  P 

(1-R)  (=i)         +  2R  cos  0         _L    -   (1+R)  =  0 

*2  L  2 


(A6) 


with  the  solution 


1   -  R 


[~-     R  cos  0    2    +7    R2  cos2012  +  (1-RZ)J  (A7) 


Dividing  (A5)  by  P-,      we  get 


M*  -3 


(A8) 
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From  (A2) 


sin 


13  ~     P. 


sin 


sin 


12 


12 


3/P: 


In  the  center-of-mass   system 


tan  6 


c.  m. 


1  2cm 


v2  sin  «u 


(A9) 


sin  0 


12 


v-,   cos   y,-,   -  V 
L  \  L  c.  m. 


cos 


c.  m. 


12 


Therefore 


Now 


c.  m. 


m 


K 


v 


i       Vp, 


m,+mT^    v-, 
a       iv      c. 


sm 


m  ,+m,.^ 
a        K 


Pl/P2 
1   +  1/R 


m 


K 


tan  Q 


12 


12cm.  P1/P2 

cos  0l2-nrT7R 


(A10) 


L^(proj)  =  L^  cos    \-3  =  L,?  cos   (  -^ —  cos        y?) 
We  find     y,     by  using  unit  vectors  and  the  three  equations 


n.     A   n7  =     p  sin 


1 


12  ' 


n,j    X    n      =  p  sin  6.  ~, 


(All) 
(A12) 

n.,  X     n?  =  p  sin  6?^  .  (A13) 

Solving  for     p     and  taking  the    j     components,    we  obtain 


(v1a2-Y2a1)  =    L_  (y^   _  Yia3) 


1 


cos  9 


cos 


(Y3a2  -   Y2a3) 

12  —-13  COS   y23 

Using  any  of  various   combinations  of  the  two  equalities,    we   solve  for 


Yl    sin  U??   -  Y?  sin  " i  -> 


Y 


sin 


(A14) 


12 
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B.      Inelastic-Scattering  Analysis  for  K~+d-»K     +  p  +  n 


1     =     incident  K 


2  =  scatter  K 

3  =  recoil  P 

A  A 

4  =  recoil  n 


Given:     P-,,    P.,,    direction  cosines  of    n.,    n?,    n. 
Find:     P.  ,  P.,    direction  cosines  of    fl.     . 


We  let     P.    =  P       i      and  therefore   P        =  P        =  0  . 

From  momentum  and  energy  conservation  we  get  the  equations 

P        =   P        +  P        +  P 

Xl  x2  X3  X4 

0         =  P        +  P        +  P 

^2  ^3  ^4 

0         =  P        +  P        +  P 


(A15) 
(A16) 
(A17) 


r>    2  p    2  p    2 

2               3                4  ^ 
+  +  +  Q, 


(A18) 


m 


K 


m 


K 


m 


m 


n 


where     Q  =  m     +m      -  m,  =   2.229. 
p  n  d 


We  immediately  can  get  P         and    P         from  Eqs.    (Al6)  and  (A17). 

y4  z4 


We  define 


K 


A       P2' 


2  P      2       P      2 

^4     +      Z4 


m 


K 


m 


m 


n 


+  2Q 


K9  =      P        +  P 
2  x2  x3 


Then,    from  Eqs.    (A15)  and  (Al6), 


Pl    =  K2  +  Px4  ■ 

Pl  X4 

=  K,    + 


m 


K 


"1  m 


n 


(A19) 
(A20) 
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Squaring  Eq.    (A19),    substituting  in  Eq.    (A20),    and  solving  the  resulting 

quadratic  we  get  ■  

/K2  +  K22-(l-^-)(K22-m     Kj) 

Px     =  -^ 2 ■    (A21) 

4  m 

(1 — ) 

m 
n 

Substituting  Eq.    (A21)  in  Eq.    (A19)  we   solve  for 


mK       __  mK     „    2       ..    mK    . 

K-,  +      / K        +  (1 )  m„  K. 

mn  2         \|  mn         z  v      mn      '       K       1 

1   "  "  

,        mK 


mn 


(A22) 
We  want  a  coordinate  system  in  which    n.     is  in  the     x    direction 
and    n-,     has  a  dip  angle  equal  to  zero.      To  get     n.     along  the     x     axis 
we  make  2  successive  transformations  .      We  rotate  about  the     z     axis 
through  an  angle     9    arid  then  about  the     y    axis  through  an  angle 

n     -  <J>): 

K      2 

n 
x 

n 

y 

n 
z 

The  given  variables  are 


f    sin  cj> 

0 

COS    (J) 

0 

1 

0 

<r  cos<t> 

0 

sin  cj> 

n         =a         n         =3,    n        =Y, 
Xl  •  ?\  Zl 

cos  <j>  =  y  s 


•  *  n — 2" 

sin  cp  =  V    1    -   y 


cos  6  -  _  , 


sin  6  = 


nI  i-y' 

p 


44. 


T 


o  get     X-,  =  0     we  rotate  about  the    x    axis  through  an  angle      £: 


vhere 


We  get 


Let 


then 


V\        /'  °  ° 

n    "      1       =10  cos   i  sin  i      I, 

\ 

n    "    /  \0  -sin  £  cos   £, 

z 


n       ■ 

Jj> A 

cos    c,  = £— =      a  , 

X-, 


■      t  Zz  A         K 

sm  £,  =  =  b 

-s/  1-n   ,2 
x-, 


n    '  =  an     +  Sn     +  v  n      , 
x  x       r    y       '     z 

n    •  =  -         r:  n    +        ■  n 

y  r      ~I  x    n        2     Y 


N/    1    -  v  \l    1 


n    '    =   --53£ n      -         ^  n      +  -/  1   -  ^  n 


v    1   -  Y  \f   1  -  y 


n        =6,n        =€,n        =   P  : 
x2  ^2  Z2 


y 


n^    '    =  a6  +  (3e    +  YP» 
(36  +  ae 


X2 


n 


Y2      /TT7 


,  _         oy6  -   (3y6    +  (1   -   Y    )P 
z 


'^ n — 7 

2  N    1   -  y 


and 


45 


n    "  -  an     +  Sn     +  an       , 
x  x       r    y  z 

(A(3  +  Bay)n     +  (Aa  -   B(3a)n 

n    "  =   - ^ ^       +  Bn 

y  2  z 

1    -   Y 

(B(3  -   Aay)n      -    (Ba  +  A(3y)n 

n    "  =    ^L ^    +  An        , 

z  ,2  z 

1    -   y 


where 


A   ■=    a.  \l  1  -  y 


B    ^  b  V  1    -  y2 


ae    - 

(36 

Vi- 

ay6 

Y       ^ 

-  Py* 

1    -  n 
> 

+  (1    - 

:2 
Y    )   P 

Jl 

2 
-   Y 

Jl 

-   n       ■ 
X2 

2 

Now  we  can  solve  for     n"  ,    and    n~"    and  then  for     n    ''     .    We  use 


the  reverse  transformation  to  get     n- 


4k 


C.      Momentum  Transfer 


In  the  center -of-mass  sytem, 


the  momentum  transfer  is 


q 


K-  kT1 


By  law  of  cosines,    we  have 


q2  =  K2  +  K'  Z  -    2  K  K'    cos   0 .  , 
n  12  c. 


m. 


For  an  elastic  scattering  we  have 

m  , 
K  =  K1    *■ — 


m  ,   +  mT. 
d  K 


PK=  0.792  PR 


Therefore  we  obtain 


o      =  2  YT  (  1   -  cos  0..,  ) 

H  x  12  c.  m.  ' 


=     nTTkn/    1    -   COS   0 


12  c.  m. 


=   1.12  P,_    4    1   -  cos  0,  , 

K  12  c.  m. 
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